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HIGHLIGHTS 


►  We  image,  segment  and  reconstruct  a  commercial  Gore  cathode  catalyst  layer. 

►  We  introduce  a  method  to  partially  or  totally  fill  the  pore  space  with  water. 

►  We  differentiate  between  a  hydrophilic  and  hydrophobic  pore  space. 

►  The  method  is  based  on  the  pore  size  distribution  of  the  pore  space. 

►  We  use  the  model  to  calculate  parameters  accessible  surface  area  and  diffusivity. 
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We  reconstruct  a  section  of  the  cathode  catalyst  layer  (CCL)  of  a  Gore  PEMFC  membrane  electrode 
assembly  three-dimensionally  with  nanometre  scale  resolution.  Subsequently,  we  present  a  new 
modelling  method  to  fill  the  pore  space  of  this  matrix  stepwise  with  water,  enabling  the  description  of 
varying  saturation  conditions  of  the  CCL.  The  method  is  based  on  a  3D  pore  size  distribution  and  enables 
to  differentiate  between  a  hydrophilic  and  a  hydrophobic  CCL.  It  utilizes  a  sequence  to  fill  the  pores 
according  to  their  size,  going  from  small  to  large  (hydrophilic)  or  vice  versa  (hydrophobic),  until  a  pre¬ 
defined  value  of  water  saturation  is  reached.  We  compare  both  cases  by  calculating  an  effective  diffu¬ 
sivity  for  oxygen  in  nitrogen  in  all  spatial  directions.  Both  the  hydrophilic  and  the  hydrophobic  case 
display  a  similar  ability  to  transport  oxygen  up  to  approximately  50%  water  saturation  of  the  pore  space. 
At  higher  water  saturation,  we  calculate  larger  diffusivity  values  for  the  hydrophobic  case.  Finally,  we 
calculate  the  specific  reaction  surface  area  that  is  accessible  from  the  gas  diffusion  layer  via  unfilled  pores 
for  all  water  saturation  conditions.  At  50%  saturation,  the  hydrophobic  case  displays  a  twenty  times 
larger  reaction  surface  area  than  the  hydrophilic  case. 

©  2012  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Optimized  and  effective  water  management  is  one  of  the  key 
challenges  that  remain  in  proton  exchange  membrane  fuel  cell 
(PEMFC)  research,  as  it  influences  the  performance  of  nearly  all 
crucial  components  [1-4].  The  ionic  conductivity  of  the  membrane 
depends  on  membrane  hydration.  On  the  one  hand,  insufficient 
water  saturation  of  the  membrane  leads  to  a  higher  resistance  and 
thus  higher  losses,  and  can  cause  total  failure  of  the  fuel  cell  due  to 
irreparable  damage  to  the  membrane.  On  the  other  hand,  excessive 
liquid  water  will  clog  both  catalyst  and  gas  diffusion  layers  (GDLs), 
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slowing  gas  diffusion  to  reaction  sites  significantly.  Consequently, 
establishing  and  maintaining  an  optimal  balance  is  mandatory. 

As  water  is  produced  in  the  cathode  catalyst  layer  (CCL),  flood¬ 
ing  is  primarily  a  problem  on  the  cathode  side  of  a  PEMFC.  While 
water  transport  mechanisms  in  the  flow  channel,  GDL  and 
membrane  have  been  well  studied  [5—9],  the  CCL  is  often 
approximated  as  an  infinitesimally  thin  layer.  In  reality,  the  CCL  is 
a  porous  medium  with  the  water-forming  reaction  occurring 
throughout  the  porous  network.  Models  that  consider  the  CCL  as 
a  truly  spatial  structure  which  is  partially  saturated  with  water  are 
still  uncommon.  Recently,  Das  et  al.  formulated  a  one-dimensional 
analytical  model  that  allows  for  the  extraction  of  liquid  water 
profiles,  describing  water  content  as  a  function  of  layer  depth,  for 
both  hydrophilic  and  hydrophobic  CCLs  [10]. 
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Nomenclature 

A 

cross-sectional  area  of  a  pore  (m2) 

C 

circumference  of  a  pore  (m) 

Dk 

effective  Bosanquet  diffusivity  (m2  s-1) 

Dbutk 

effective  bulk  diffusivity  (m2  s_1) 

Df 

effective  Knudsen  diffusivity  (m2  s  j) 

jj  bulk,0 

free,  substance-specific  bulk  diffusivity  (m2  s-1) 

jjl<n,0 

free,  substance-specific  Knudsen  diffusivity  (m2  s_1) 

nbulk  ,* 
uk 

dimensionless  bulk  diffusion 

DKnf 

dimensionless  Knudsen  diffusion 

I< 

curvature  (m-1) 

I<n 

Knudsen  number 

l 

physical  length  scale  of  the  porous  medium 
involved  (m) 

Pc  ap 

capillary  pressure  (N  m-2) 

Vmean 

thermal  velocity  (m  s-1) 

6 

contact  angle  (°) 

K 

subscript  defining  spatial  direction 

A 

mean  free  path  (m) 

a 

surface  tension  (N  m-1) 

On  the  one  hand,  water  in  the  CCL  is  necessary  to  assist  proton 
conduction  to  reaction  sites  and  is  thus  beneficial  for  optimal 
operation.  On  the  other  hand,  water-flooded  pores  will  significantly 
hinder  oxygen  diffusion  to  reaction  sites  and  thus  have  a  negative 
influence  on  efficiency.  In  other  words,  one  can  assume  that  there  is 
an  optimum  water  content  for  the  CCL. 

Fundamentally,  the  liquid  water  content  of  the  CCL  is  governed 
by  the  balance  of  evaporation  and  condensation  of  water  as  well  as 
the  pressure  distribution  in  the  porous  medium.  In  the  following 
we  only  focus  on  the  pressure  distribution  which  can  be  described 
by  the  capillary  pressure  in  the  pore  network  and  the  pressure  on 
the  interfaces  between  the  CCL  and  membrane,  and  the  CCL  and 
GDL.  The  capillary  pressure  pc ap  at  an  interface  between  a  gas  and 
a  liquid  is  determined  by  the  Young-Laplace  equation: 

Pc  ap  =  &K  (1) 

with  the  surface  tension  a  between  gas  and  liquid  and  the  curva¬ 
ture  I<  of  the  interface.  While  the  surface  tension  -  assuming  an 
approximately  homogeneous  temperature  and  material  distribu¬ 
tion  throughout  the  CCL  —  is  a  constant,  the  curvature  varies 
strongly  depending  on  the  actual  pore  dimensions  at  the  interface 
location  and  the  contact  angle  at  the  associated  triple-phase  line. 
Additionally,  a  porous  media  usually  is  not  composed  of  a  bundle  of 
unconnected  capillaries  but  instead  consists  of  a  complex  network 
of  interconnected  pores.  Accordingly,  equilibrium  pressure  must  be 
established  at  each  intersection,  which  is  influenced  by  and  itself 
influences  all  connected  pore  space.  Considering  this  degree  of 
complexity,  it  is  hardly  surprising  that  until  now,  only  methods  that 
somewhat  simplify  the  structure  before  offering  a  solution  have 
been  reported.  There  are  many  different  approaches  to  accomplish 
this.  In  the  porous  media  utilised  in  fuel  cells,  these  models  have 
mainly  been  applied  to  the  GDL  to  date.  This  includes  models  that 
generalize  the  pore  structure  (e.g.  pore  network  model  [11,12])  or 
describe  the  layer  with  parameters  such  as  effective  porosity  and 
permeability  [13,14]. 

In  recent  years,  scanning  electron  microscopy  tomography 
(SEMt)  [15]  was  utilized  to  study  the  microstructure  of  solid  oxide 
fuel  cells  (SOFCs).  Electrode  material  was  removed  layer  by  layer  by 
a  focused  ion  beam  (FIB)  and  images  of  each  layer  were  combined 


to  achieve  a  three-dimensional  reconstruction  of  the  electrodes 
investigated  [16-20].  This  method  was  also  applied  to  lithium-ion 
batteries  [21-23]  and  PEMFC  CCLs  [24,25].  In  PEMFC  CCLs,  this 
allowed  for  a  detailed  differentiation  of  pore  and  bulk  material. 

In  the  reported  work,  we  combine  the  3D  SEMt  reconstruction 
of  a  CCL  with  stepwise  filling  of  the  pores  with  water  for  both 
a  hydrophilic  and  a  hydrophobic  CCL.  The  method  is  explained  and 
validated.  We  investigate  how  the  parameters,  diffusivity  and 
accessible  surface  area,  change  as  the  water  content  increases. 

2.  Experimental 

The  cathode  of  a  commercially  available  Gore  PRIMEA  A510.1 
M710.18C510.4  PEMFC  membrane  electrode  assembly  (MEA)  was 
studied  with  SEMt.  Spatial  directions  are  defined  as  follows  and 
maintained  in  all  further  discussion:  The  z  direction  is  the  general 
translation  direction  of  the  FIB  whereas  the  x  and  y  directions 
define  the  planes  visible  in  all  SEMt  images.  The  y  direction  indi¬ 
cates  the  through-plane  direction  from  the  interface  between  the 
CCL  and  membrane  to  the  interface  between  the  GDL  and  CCL 
(Fig.  1). 

A  Zeiss  ‘Neon  40EsB’  was  used  for  the  SEMt.  A  2.5  x  2.5  pm2 
platinum  protection  layer  was  deposited  on  top  of  the  porous 
cathode  catalyst  layer  with  a  30  kV  acceleration  voltage  (AV)  and 
a  current  of  5  pA  for  200  s.  Two  reference  lines,  one  orthogonal  and 
the  other  with  an  angle  of  55°  relative  to  the  cutting  plane,  were 
inscribed  into  the  platinum  layer  with  the  FIB  at  30  kV  AV  and 
50  pA,  providing  a  method  to  determine  slice  thickness  which  is 
independent  of  surface  skew  or  irregularities.  A  cavity,  a  rectan¬ 
gular  prism  of  dimensions  6  pm  x  4  pm  x  10  pm  (x,  y  and  z 
directions),  was  removed  from  the  sample  by  FIB  milling  at 
30  kV  AV  and  500  pA.  As  a  starting  point,  the  sidewall  of  the  cavity 
in  the  z  direction  was  finely  polished  at  30  kV  AV  and  50  pA.  Then 
the  FIB  removed  a  thin  layer  of  the  sample  in  the  z  direction  at 
30  kV  AV  and  50  pA,  while  the  SEM  generated  an  image  of  every 
slice  between  ablation  steps  at  2  kV  AV.  The  process  was  subse¬ 
quently  repeated  to  yield  a  series  of  120  images  with  an  average 
spacing  of  13.6  nm  between  images  and  a  resolution  of  2.5  nm  per 
pixel  (Fig.  1 ).  The  images  were  carefully  assessed  for  damage  to  the 
cutting  plane  by  either  the  slicing  or  imaging  step.  No  indication  of 
deformation  on  the  cutting  plane  could  be  detected.  The  pore  walls 
behind  the  cutting  plane,  visible  due  to  the  depth  of  focus,  dis¬ 
played  no  change  in  morphology  between  successive  images. 
While  a  small  degree  of  distortion  can’t  be  ruled  out  due  to  FIB/SEM 
being  an  invasive  method,  this  error  is  thus  judged  to  be  negligible. 

The  complete  experimental  reconstruction  procedure  of  the 
SEMt  including  alignment  was  performed  according  to  the 
approach  described  in  [21,25].  The  segmentation  was  accomplished 
semi-automatically  by  first  applying  a  threshold  and  then  manually 
correcting  the  images. 

Methods  such  as  calculating  the  pore-size  distribution  according 
to  the  approach  employed  in  this  work  demand  cubic  voxels. 
Accordingly,  the  resolution  of  the  x  and  y  direction  was  first  rescaled 
to  3.4  nm  x  3.4  nm  by  linear  interpolation,  resulting  in  a  voxel  size  of 
3.4  nm  x  3.4  nm  x  13.6  nm.  Then  each  voxel  was  divided  into  4  sub¬ 
voxels  measuring  3.4  nm  x  3.4  nm  x  3.4  nm.  The  final  geometry  was 
cropped  to  a  cube  of  1.632  pm  x  1.632  pm  x  1.632  pm  and 
480  x  480  x  480  voxels.  A  preliminary  porosity  of  58.08%  was 
calculated. 

On  every  image,  small  non-connected  areas  attributed  to  errors 
during  the  semi-automatic  segmentation  process  were  eliminated 
from  both  the  solid  and  porous  phases.  This  was  performed  by 
applying  the  bwareaopen  function  in  Matlab  [26]  in  2D  which 
removed  all  non-connected  2D  particles  with  a  size  smaller  than  10 
pixels  set  to  an  8-connected  neighbourhood  (all  pixels  that  share  an 
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Fig.  1.  a)  Schematic  diagram  of  the  experimental  set-up.  A  small  sample  of  a  PEMFC  cathode  is  ablated  slice  by  slice  with  a  focused  ion  beam.  A  scanning  electron  microscope  angled 
at  36°  relative  to  the  sample  surface  plane  creates  one  image  per  slice.  Schematic  as  published  in  [21].  b)  One  of  the  images  recorded  between  ablation  steps.  The  cutting  plane  and 
sample  surface  with  two  reference  lines  imprinted  into  a  deposited  platinum  layer  for  the  alignment  process  are  visible.  Due  to  the  camera  angle,  images  must  be  aligned. 


edge  or  vertex  are  defined  as  connected).  The  result  is  an  inter¬ 
mediate  porosity  of  58.12%. 

The  bwconncomp  function  in  Matlab  [26]  was  used  in  3D  to 
remove  small  non-connected  3D  areas  and  evaluate  differences  in 
3D  connectivity.  This  was  performed  by  the  following  steps  for  a  6- 
connected  (voxels  must  share  a  face  to  be  connected),  18-connected 
(voxels  must  either  share  a  face  or  edge  to  be  connected)  and  26- 
connected  (voxels  must  either  share  a  face,  edge  or  vertex  to  be 
connected)  neighbourhood  respectively: 

i)  The  largest  connected  solid  region  in  the  3D  geometry  is 
identified.  Depending  on  connectivity  definition,  this  region 
consists  of  98.78%  (6-connected),  98.80%  (18-connected)  and 
98.80%  (26-connected)  of  the  voxels  defined  as  solid  in  the 
geometry. 

ii)  All  other  solid  regions  not  connected  to  this  largest  region  are 
converted  to  pore  space. 

iii)  The  largest  connected  pore  in  the  3D  geometry  is  identified. 
Depending  on  connectivity  definition,  this  pore  consists  of 
between  99.88%  (6-connected),  99.90%  (18-connected)  and 
99.91%  (26-connected)  of  the  voxels  defined  as  pore  in  the 
geometry. 

iv)  All  other  pores  not  connected  to  this  largest  pore  are  con¬ 
verted  to  solid. 

The  results  are  porosities  of  58.14%  (6-connected),  58.15%  (18- 
connected)  and  58.15%  (26-connected).  The  6-connected  porous 
3D  matrix  was  chosen  as  starting  point  for  all  further  investigations 
(Fig.  2). 

3.  Mechanism  of  water  filling 

One  of  the  most  important  parameters  for  characterising 
a  partially  flooded  porous  medium  is  the  wettability  of  the  pore 
walls.  Due  to  the  small  dimensions  of  the  CCL,  measuring  an 
internal  contact  angle  is  not  yet  feasible.  Most  authors  characterise 
the  catalyst  layer  by  placing  a  static  droplet  on  the  surface  of  a  CCL 
to  gain  at  least  a  general  indication  concerning  the  wettability  [27— 
31].  The  disadvantage  of  this  approximation  is  that  there  is  no 
guarantee  that  the  generated  values  are  applicable  to  the  pore 
space  consisting  of  pores  on  a  totally  different  scale.  Additionally, 
a  certain  amount  of  contact  angle  hysteresis  is  to  be  expected, 
further  increasing  the  possible  error.  To  reduce  this  uncertainty,  we 
consider  the  equation  first  introduced  by  Finn  [32]: 

J  Pap  =  J  ffCOS  6  (2) 

A  C 


sectional  area  A  without  knowledge  of  the  exact  interface  shape. 
At  6  =  90°,  the  capillary  pressure  changes  sign,  thus  altering  the 
fundamental  behaviour  concerning  where  liquid  water  will  pref¬ 
erably  reside  to  minimize  pressure.  In  other  words,  in  the  hydro¬ 
philic  case,  small  pores  will  primarily  fill  with  water;  in  the 
hydrophobic  case,  the  larger  pores  will  be  more  attractive  to  water. 
Our  model  describes  these  two  different  cases. 

A  three-dimensional  pore  size  distribution  (PSD)  naturally 
contains  information  that  differentiates  between  different  sizes  of 
pores.  Consequently  we  model  the  fundamental  behaviour  of 
hydrophilic  and  hydrophobic  water  filling  by  first  calculating 
a  spatially  resolved  PSD  and  then  gradually  filling  the  pore  space 
with  water.  This  water  filling  depends  solely  on  pore  size.  In  the 
hydrophilic  case  the  smallest  pores  are  filled  first  and  then  larger 
and  larger  pores  are  filled  until  a  predefined  threshold  of  water 
saturation  is  reached  (Fig.  3).  This  was  performed  in  10%  steps  of 
total  pore  volume,  until  all  pore  space  is  filled.  In  the  hydrophobic 
case,  the  mechanism  is  inverted.  The  largest  pores  are  filled  first 
and  then  smaller  and  smaller  pores  are  filled  up  to  a  predefined 
threshold.  This  was  also  performed  in  10%  steps  of  pore  space  total 
pore  volume,  until  all  pore  space  is  filled.  To  accomplish  this,  we  use 
an  adaptation  of  the  approach  which  was  described  in  reference 
[33]  and  further  developed  in  reference  [25].  In  a  first  step,  every 
voxel  belonging  to  the  pore  space  of  the  3D  reconstruction  matrix  is 
identified.  Then  the  distance  to  the  nearest  pore  wall  is  calculated. 
Spheres  with  a  radius  of  this  minimum  distance  are  inscribed  into 
each  voxel  of  the  pore  space,  with  larger  spheres  overwriting 


,  ,  .  1 1  .  ^  r  ,  |  .  ,  Fig. 2-  Visualization  of  the  3D  reconstruction  of  a  PEMFC  cathode  used  as  starting 

where  the  capillary  pressure  Pcap  over  an  interface  can  be  calculated  point  for  water  modelling  in  the  reported  work.  The  cubic  representation  of  the  bulk 

with  the  contact  angle  6  along  the  circumference  C  and  the  cross-  phase  measures  1.632  pm  along  each  edge. 
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a 


empty  channel 


b 


pore  size  distribution 


c 


hydrophilic  (50%  water) 


hydrophobic  (50%  water) 


Fig.  3.  a)  Empty  channel  (white  =  pore,  black  =  bulk),  b)  Spheres  inscribed  into  the 
pore  as  a  measure  for  the  pore  size  distribution  (PSD).  Smaller  spheres  are  depicted  in 
darker  grey.  Larger  spheres  overwrite  smaller  spheres,  c)  Hydrophilic  pore  with 
a  water  saturation  of  50%  (white  =  pore,  black  =  bulk,  grey  =  water).  Small  spheres  of 
the  PSD  are  filled  first,  d)  Hydrophobic  pore  with  a  water  saturation  of  50%  (white- 
=  pore,  black  =  bulk,  grey  =  water).  Large  spheres  of  the  PSD  are  filled  first. 

smaller  spheres.  This  avoids  unphysical  gaps  between  spheres,  as 
would  occur  in  a  packing  of  spheres  without  overlap.  A  grey-scale 
value  of  100  describing  water  (with  0  -  black-and  255  -  white¬ 
depicting  pore  and  bulk  material  respectively)  is  assigned  to  this 
resulting  set  of  spheres.  The  geometry  is  filled  in  10%  steps  of  pore 
volume  from  smallest  to  largest  (hydrophilic  case)  or  largest  to 
smallest  (hydrophobic  case)  (Fig.  4,  Fig.  5).  To  verify  the  method, 
a  PSD  [25,33]  of  the  remaining  pore  space  is  calculated  with  Geo- 
Dict  [34]  and  depicted  in  Fig.  6.  The  stepwise  water  flooding  was 
implemented  in  Matlab  [26]. 

This  three-dimensional  PSD  water  model  does  not  consider 
water  transport  dynamics,  such  as  contact  angle  hysteresis,  or  the 
possibility  of  capillary  stops  due  to  a  varying  pore  radius  along 
a  capillary.  However,  we  describe  an  important  base  case  where  the 
surface  energy  of  the  geometry  is  at  its  minimum  for  a  specific 


degree  of  water  saturation.  This  procedure  lets  us  take  full  advan¬ 
tage  of  the  3D  reconstruction  of  the  real  microstructure  presented 
in  this  work  for  the  calculations  of  parameters  in  dependence  on 
water  content,  without  having  to  simplify  pore  structure. 

4.  Characterization  of  the  partially  filled  3D  reconstruction 

The  diffusion  rate  of  oxygen  in  air  is  approximately  105  times 
higher  than  in  water  [35].  As  a  result,  we  assume  that  only  the  pores 
that  are  not  filled  with  water  contribute  to  the  oxygen  supply  of  the 
catalyst  layer.  This  allows  us  to  graphically  add  the  water  phase  to  the 
bulk  phase  for  every  degree  of  water  saturation  of  our  3D  recon¬ 
struction,  with  the  colour  black  describing  pore  space  and  the  colour 
white  depicting  the  sum  of  water  and  bulk  phase,  thus  significantly 
reducing  calculation  time  for  the  following  calculations. 

4.1.  Calculation  of  gas  diffusivity  values  in  dependence  on  water 
saturation 

The  appropriate  diffusion  mechanism  in  a  porous  medium  is 
determined  by  the  Knudsen  number  Kn : 

I<n  =  J  (3) 

X  is  defined  as  the  mean  free  path  between  two  successive  inter- 
molecular  collision  events  and  is  a  specific  value  for  a  given  gas  at 
a  constant  pressure  and  temperature.  I  describes  a  representative 
geometric  physical  length  scale  of  the  porous  medium  involved, 
which  is  typically,  but  not  necessarily,  an  average  pore  diameter.  If 
this  physical  length  scale  is  much  smaller  than  A,  Kn  will  become 
very  large  and  instead  of  intermolecular  collisions,  molecules  will 
predominantly  collide  with  pore  walls.  This  type  of  diffusion  is 
called  Knudsen  diffusion  and  can  be  modelled  by  a  random  walk 
approach  [36-38]. 

For  small  Knudsen  numbers,  diffusion  mainly  occurs  due  to 
collisions  between  molecules.  This  type  of  diffusion  is  called  bulk 
diffusion  and  is  typically  modelled  by  the  well-known  continuum 
mechanics  approach. 

At  intermediate  Knudsen  numbers,  as  are  applicable  to  the 
reconstructed  CCL,  both  intermolecular  and  molecule-wall  colli¬ 
sions  contribute  significantly.  In  this  case,  the  overall  effective 
diffusivity  can  be  approximated  by  Bosanquet’s  formula  [37,39]: 


Fig.  4.  2D  sectional  view  of  0.2  jam  edge  length  in  z  direction  for  varying  filling  levels,  a)  0%,  20%  and  60%  water  content  of  a  hydrophilic  3D  reconstruction,  b)  0%,  20%  and  60%  water 
content  of  a  hydrophobic  3D  reconstruction. 


264 


T.  Hutzenlaub  et  al.  /  Journal  of  Power  Sources  227  (2013)  260-266 


Fig.  5.  Visualization  of  a  partially  filled  3D  reconstruction.  20%  of  the  pore  space  is  filled  with  water  (blue),  a)  Hydrophilic  case,  small  pores  are  filled  first,  b)  Hydrophobic  case,  large 
pores  are  filled  first.  (For  interpretation  of  the  references  to  colour  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  this  article.) 
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Fig.  6.  Pore  size  distribution  (PSD)  of  the  remaining  pore  space  at  varying  filling  levels, 
a)  Hydrophilic  filling;  small  pores  are  filled  first.  Accordingly,  the  PSD  distribution  is 
cropped  stepwise  from  left  to  right.  Irregularities  are  due  to  edge  effects  at  the 
interface  between  water,  solid  and  gas,  thus  creating  additional  possibilities  to  inscribe 
small  spheres  due  to  the  curvature  near  the  triple-phase  line,  b)  Hydrophobic  filling; 
large  pores  are  filled  first.  Accordingly,  the  PSD  distribution  is  cropped  stepwise  from 
right  to  left.  Irregularities  are  due  to  edge  effects  at  the  interface  between  water,  solid 
and  gas,  thus  creating  additional  possibilities  to  inscribe  small  spheres  due  to  the 
curvature  near  the  triple-phase  line. 


Dk  =  ((Df )  1  k  =  x,y,z  (4) 

Both  the  effective  bulk  diffusivity  D^ulk  and  the  effective  Knudsen 
diffusivity  D%n  can  be  understood  as  the  result  of  dimensionless 
diffusivity  matrices  Dbulk*  and  £>£"*,  depending  only  on  the 
geometry  of  the  porous  media,  multiplied  by  free,  substance- 
specific  diffusivity  values  Dbulk,°  and  DKn,°,  respectively. 

Dbulk  =  Dbulk,*Dbulk,°  k  =  x,y,z  (5) 

D/Cn  =  jjKn*  DKn,0  K  =  xyz  (6) 

Dbulk’*  and  D £n’*  are  calculated  as  described  by  Becker  et  al.  [37] 
with  GeoDict  [34],  where  Dbulk;*  is  the  result  of  solving  Fields  law 
and  D^n’*  is  derived  by  utilizing  a  random  walk  approach. 
Dbulk  *D^n*  and  the  characteristic  length  l  which  is  the  average 
distance  between  molecule-wall  collisions  and  is  the  result  of  the 
random  walk  simulation  are  depicted  in  Figs.  7-9.  Dbulk,°  is  defined 
as  either  the  self-diffusion  coefficient  of  a  pure  gas  or-in  the  case  of 
a  binary  gas  mixture  such  as  oxygen  in  nitrogen  as  we  can  assume 
in  the  CCL-the  measurable  diffusion  coefficient  in  free  space.  At 
standard  conditions  (25  °C,  191.3  kPa)  Dbulk-°  =  20.86  mm2  s_1  [37]. 
DKn,°  is  given  by: 

DKn’°  =  Ifo  mean  (7) 

At  standard  conditions  the  mean  thermal  velocity  is 
Vmean  =  444.1  m_1  [37].  Thus  we  are  now  able  to  determine  the 
effective  diffusivity  coefficient  DK  of  oxygen  in  nitrogen  in  the  CCL  in 
dependence  on  the  water  content  of  the  pores  with  Bosanquet’s 
formula.  The  result  is  depicted  in  Fig.  10. 

Both  the  hydrophilic  and  the  hydrophobic  case  display  a  similar 
ability  to  transport  oxygen  up  to  approximately  50%  water  satura¬ 
tion  of  the  pore  space.  At  higher  water  saturation,  we  calculate 
larger  diffusivity  values  for  the  hydrophobic  case.  The  pore  space  in 
the  hydrophilic  case  is  entirely  blocked  at  approximately  70%  water 
saturation.  In  the  hydrophobic  case,  connectivity  through  the 
sample  is  provided  up  to  a  value  of  approximately  90%  water 
saturation. 
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Fig.  7.  Dimensionless  bulk  diffusivity  in  all  spatial  directions  for  both  the  hydrophilic 
and  the  hydrophobic  cases.  Connecting  lines  between  data  points  are  added  for  better 
visibility. 
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Fig.  8.  Dimensionless  Knudsen  diffusivity  in  all  spatial  directions  for  both  the  hydro¬ 
philic  and  the  hydrophobic  cases.  Connecting  lines  between  data  points  are  added  for 
better  visibility. 

4.2.  Calculation  of  specific  GDL-connected  surface  area  in 
dependence  on  water  saturation 

In  general,  the  effective  diffusivity  in  all  spatial  directions  gives 
a  good  indication  of  how  the  reconstructed  part  of  the  CCL 


Fig.  10.  Effective  diffusivity  in  all  spatial  directions  for  both  the  hydrophilic  and  the 
hydrophobic  cases.  In  the  hydrophilic  case,  diffusion  through  the  reconstructed  section 
comes  to  a  complete  standstill  at  approximately  70%  water  content  in  the  pore  space. 
In  the  hydrophobic  case,  diffusion  through  the  reconstructed  section  is  possible  until 
approximately  90%  water  content.  For  both  cases  all  three  spatial  directions  yield 
similar  results.  Connecting  lines  between  data  points  are  added  for  better  visibility. 


contributes  to  blockage  or  passage  of  oxygen  in  dependence  on 
water  saturation,  thus  characterising  its  degree  of  connectivity. 
However,  this  is  not  an  exact  measure  of  how  the  CCL-represented 
by  the  reconstructed  3D  cube-losses  reactivity  as  the  water 
content  rises.  To  quantify  this,  one  can  assume  that  the  most 
important  path  for  oxygen  diffusion  is  from  the  GDL  side  into  the 
CCL.  Thus,  we  fill  all  pores  that  are  not  connected  to  this  side  of  the 
reconstructed  3D  geometry  with  bulk  material  by  a  tool  based  on 
the  bwconncomp  function  in  Matlab  [26]  and  then  calculate  the 
remaining  surface  area  that  is  not  covered  by  water.  The  surface 
area  calculation  was  performed  by  employing  a  method  first 
introduced  by  Ohser  et  al.  [40]  with  GeoDict  [34]  (Fig.  11)  for  both 
the  hydrophilic  and  the  hydrophobic  cases  at  varying  values  of 
water  saturation.  This  surface  area  is  related  to  the  volume  of  the 
3D  reconstruction  to  determine  a  specific  surface  area  that  is 
available  for  reaction.  The  hydrophobic  case  demonstrates  consid¬ 
erably  higher  tolerance  to  water  saturation.  At  50%  saturation,  we 
calculate  a  specific  reaction  surface  area  that  is  approximately 
twenty  times  higher  than  for  the  hydrophilic  case. 
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Fig.  9.  The  characteristic  length  for  both  the  hydrophilic  and  the  hydrophobic  cases.  It 
is  defined  as  the  average  distance  between  molecule-wall  collisions  and  is  the  result  of 
the  random  walk  simulation.  Connecting  lines  between  data  points  are  added  for 
better  visibility. 


Fig.  11.  Specific  surface  area  connected  to  the  gas  diffusion  layer  side  of  the  recon¬ 
structed  sample  via  pore  space  that  is  not  filled  with  water.  At  50%  saturation,  the 
specific  reaction  surface  area  is  approximately  twenty  times  higher  than  for  the 
hydrophilic  case.  Connecting  lines  between  data  points  are  added  for  better  visibility. 
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5.  Conclusion  and  outlook 

Firstly,  we  successfully  demonstrated  a  mechanism  to  partially 
fill  a  3D  reconstruction  of  a  porous  medium  with  water,  without 
needing  to  streamline  the  geometry  to  a  more  calculation-friendly 
configuration  such  as  packed  spheres.  This  is  valuable  to  achieve 
a  more  realistic  starting  point  for  porous  media  simulations. 

Secondly,  we  compared  a  hydrophilic  and  a  hydrophobic  CCL  by 
calculating  an  effective  diffusivity  of  oxygen  in  nitrogen  and 
a  specific,  GDL-connected  surface  area.  While  diffusivity  values  of 
both  cases  are  similar,  the  surface  area  strongly  favours  a  hydro- 
phobic  CCL  as  the  layer  that  is  expected  to  produce  higher  power 
densities,  due  to  better  accessibility  of  the  catalyst  for  oxygen. 

Comparison  of  the  proton  conductivity  for  the  two  configura¬ 
tions  is  still  an  object  for  future  research. 

The  water-filled  3D  reconstruction  presented  in  this  publication 
can  be  provided  on  request. 
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